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Abstract 

Using the semiclassical approach to the energy levels of the fermions 
bound to the vortex core, we found the difference between the states 
in nonchiral and chiral superconductors, determined by the Berry 
phase. The bound states of fermions in the singly quantized vortex in 
the layered superconductor with the symmetry of superfluid ^He-A is 
E = nujo (as distinct from E = (n + l/2)a;o in s-wave superconductor) 
and thus contains the state with exactly zero energy. This is in accor- 
dance with the result obtained in microscopic theory Using this 
approach for calculations of the effect of impurities on the spectrum 
of bound state we reproduced the Larkin-Ovchinnikov result for single 
impurity in s-wave vortex the spectrum has the double period 

AE = 2ujQ and consists of two equdistant sets of levels. The same 
approach however shows that the single impurity does not change the 
spectrum E = nu)Q in the p-wave vortex . 



1 Introduction 



The low-energy fermions bound to the vortex core play the main role in the 
thermodynamics and dynamics of the vortex state in superconductors and 
Fermi-superfluids. The spectrum of the low-energy bound states in the core 
of the axisymmetric vortex with winding number = ±1 in the isotropic 
model of s-wave superconductor was obtained by Caroli, de Gennes and 
Matricon Q: 



This spectrum is two-fold degenerate due to spin degrees of freedom. The 
integral quantum number n = is a modified angular momentum of the 
bound state fermions, while the level spacing ujo is the magnitude of the 
angular velocity of the fermions orbiting about the vortex axis. The direction 
of rotation is determined by the sign of the winding number N of the vortex. 

The level spacing is small compared to the energy gap of the quasiparticles 
outside the core, uo <C A. So, in many physical cases the descreteness of 
n can be neglected. In such cases the spectrum crosses zero energy as a 
function of continuous angular momentum L^. So, one has the fermion zero 
modes. The fermions in this ID "Fermi liquid" are chiral: the positive energy 
fermions have a definite sign of the angular momentum n. 

The motion of the vortex perturbs the fermions in the vortex core and 
leads to the spectral flow of the fermionic levels from the negative energy 
vacuum to the positive energy world of excitations forming the heat bath, or 
the normal component. This spectral flow leads to the momentum exchange 
between the vortex core and the heat bath and thus to the additional force 
acting on the vortex, when it moves with respect to the heat bath |^, ^, |^, |]. 
This force was calculated in the microscopic theory in the earlier papers by 
Kopnin and coathors 0. 

Later the other types of vortices have been found to exist in superfiuid 
■^He and in high-temperature superconductors: with different winding num- 
bers and in the nonsymmetric environment and/or with the spontaneously 
broken symmetry in the core. The microscopic description of the bound-state 
fermions in such distorted vortices becomes complicated. So one needs in the 
simple phenomenological theory, which describes the low-energy motion of 
the quasiparticles in the vortex core in the same manner as the low-energy 
fermions are described in the Landau theory of the Fermi liquid. 
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Such theory was constructed in the papers |T0|, |TT| . In this approach the 
fast radial motion of the fermions in the asymmetric vortex core is integrated 
out and one obtains only the slow motion corresponding to the fermions 
zero modes. The slow low-energy dynamics is described quasiclassically in 
terms of the pair of canonically conjugated variables: the angle 9, which 
specifies the direction of the linear momentum of the propagating fermion, 
and impact parameter (or angular momentum L^). In this description the 
fermion zero modes are represented by one or several bracnches E^{Lz, 9) = 
—uJs_{9){Lz — Lzs_{9)) of the fermionic spectrum. They cross zero energy as a 
function of the continuous angular momentum L^. The number of fermion 
zero modes does not depend on 9 and is determined by the vortex winding 
number iV [0: it is 2N if the spin degrees of fredom are taken into account. 
The system of these chiral fermions form the set of the one dimensional 
Fermi-liquids. 

The final step in the calculation of the quantum spectrum of the fermions 
in the asymmetric core is the quantization of the remaining slow motion. 
It is obtained by the Bohr-Sommerfeld rule for the canonically conjugated 
variables [1^, 0. It appears that even for the asymmetric vortex the main 
structure of the quantum energy levels remains robust. One obtains one or 
several branches of the discrete energy levels. On each branch the levels 
are equidistant as in the simplest case of the axisymmetric vortex in s-wave 
superconductor. 

Here we show that the Berry phase is important for quantization. Due to 
the Berry phase the bound states of fermions in the singly quantized vortex 
in the chiral layered superconductor with the symmetry of superfiuid '^He-A 



(possible candidate is superconducting Sr2Ru04 |jT2|) is £^ = tt-cuq, while the 



s-wave superconductor one has E = {n + 1/2)ujq. Thus the spectrum of the 
p-wave vortex contains the state with exactly zero energy, which confirms the 
microscopic calculations in [|[]. We also consider the effect of impurities on 
the spectrum of bound state. We reproduce the Larkin-Ovchinnikov result 
for single impurity in s-wave vortex p|, R[ and extend it to the p-wave case. 
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2 Boundary conditions for quantization of bound 
states in the core of s-wave vortex. 



Let us start with the A^-quantum nonaxisymmetric vortex in conventional 
superconductors, where the asymmetry is caused, say, by crystal potential or 
by external perturbation. In the semiclassical approximation the quasipar- 
ticle states are characterized by their trajectories. Trajectories are straight 
lines along the direction of the linear momentum in transverse plane q^, 
with the total momentum being on the Fermi surface, q^ + qu. = p'f- Let us 
consider somewhat more general case of distorted vortex in s-wave supercon- 
ductor. The order parameter is chosen as A(r, 0) = |A(r, 0)|e***^''^''. Here 
is azimuthal angle; $(0) is the phase of the order parameter, which winding 
around the vortex is 

where $(0) is a periodic function of (f). The Bogoliubov Hamiltonian for 
quasiparticles in the vocinity of such a vortex is 

n = p- vr3 + ri|A(r,0)|cos<l>(0) + r2|A(r,0)|sin$(0) , v = — . (3) 

m 

Let us first neglect the conventional scattering on vortices, leaving only 
the Andreev reflection. In this case, in the quasiclassical limit there is no 
transition between different trajectories: quasiparticle is moving along the 
same trajectory, i.e. with the same q_L, and experience the Andreev scattering 
which does not change the momentum, but chnges the sign of the group 
velocity. Each trajectory is characterized by the angle 6 in the transverse 
plane and by the impact parameter b: 

b = r sin , q_L = x cos + y sin 6 . (4) 

Let us make two transformations of the Eq. (Q) for given q^ . The first one 
is the coordinate transformation which orients the axis x along the trajectory 

n = ~ivT3^s + Tl\A{r,(f)-e)\cos^(f)-e)+T2\A{r,(j)-e)\sm^(f)-e) . (5) 
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In a new coordinate system px = —ids, where s = rcos0 is the coordinate 
along the trajectory. 

The second transformation is intended to make the phase of the order 
parameter the periodic function of 6: 

^ ^ ^iNr,e/2^^-iNr,e/2 ^ (7) 

-iwTs^, + ri|A(r,0)|cos(l'(0) + iV0) +r2|A(r,0)|sin(l>(0)+Ar0) , (8) 

= 0-^. (9) 

Note that for the axisymmetric vortex there is no more dependence on 6: 

n = -ivT-ids + Ti\A{r)\ cos {N(l)) + T2\A{r)\ sin {N(l)) . (10) 

The quantity N6/2 plays the part of the Berry phase: when one changes the 
direction of the trajectory, the only change in the wave function is multipli- 
cation by e*^"'3^/2. In more general case of nonsymmetric vortex this Berru 
phase is the adiabatic change in the phase of the wave function, when one 
continuously rotates the trajectory of the quasiparticle by the angle 9 in the 
transverse plane. 

Quantization of the motion along s gives rise to generally fermion zero 
modes of the type: 

E{L,,e) = -LUoie){L,~L^^\e)) (ii) 

The next step is quantization of the azimuthal motion. Here it is important 
that the wave function in terms of 6 has the boundary condition 

^(0 + 27r) = (-l)^^(^) (12) 

which follows from Eq.(^. 

The shift L^°^9) from zero is antisymmetric because of the "CPT"- 
theorem: 

L(°)(^) = -Lf(0 + 7r) (13) 

Below Eqs(|Tl],|T3|) will be explicitly shown for a weakly asymmetric vortex. 

The quantum Hamiltonian for the canonically conjugated variables Lz = 
—ide and 9 is 

n = -^MO), (-^-4°^^^))} , (14) 
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where { , } is anticommutator. The normahzed wave function of the eigen 
state with the energy E is 



xijm) = (—--) exp U / de' i—— + L^'>\9')\\ . (15) 



The energy eigenvalues 7Y\I/„(^) = i?„\&„(6') are found from the requirement 
that according to Eg . ([121) the wave function changes sign after encir- 

chng the origin in the momentum space 6 for odd N and does not change if 
is even, i.e. the phase of the wave function changes by 7r(2n + A^). This 
gives 

/''^ Ji- + dOL^^^e) = -27r (n + -] . (16) 



JO uJo{0) Jo V 2 

The second term in the Ihs of Eq. ([T6| ) is zero according to Eq. (^) , and one 
obtains the quantization which is essentially the same as for the axisymmetric 
vortex: 



,^0 (^), 



ft 

En = - , , even N. (18) 



For even there is a state with zero energy. 



3 Fermion zero modes in a weakly asymmet- 
ric vortex with odd N in 5-wave supercon- 
ductors. 

Let us consider a weakly asymmetric vortex in which the anisotropy is small, 
say, |$(0)| -C 1. In the extreme limit, when 6 = and $(0) = 0, the 
trajectory with = for s > and = vr for s < gives exactly zero energy 
if is odd. This is becuase the order parameter along this trajectory is real 
and changes sign when the origin is crossed. The index theorem requires that 
there is the normalizable eigen state with zero energy. 
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Let us choose an odd A^, then for small b the perturbation theory can be 
constructed. Since (j) is close to or tt one has 

= 0- ^^1 (1- sign s)-9 . (19) 

and 

n = n'-^^ + , (20) 

= -tvTs-^^+n\A{s,e)\ sign s , (21) 

A(s,^)=A(^|s|,^(l-signs)-^) (22) 

7^(1) = r2|A(s,^)|sin(<l(0) + A^0) ^ (23) 

^ r2|A(s,^)| I^AtA + (sign s) 4' (| (1 " sign «) - • (24) 

The Hamiltonian TC^^^ has fermion zero mode, \l/*^°^(s), so that the energy is 
the average of Ti.^^^ over this mode: 

E{b,e) = (t^w 

-iV6(^M^^^±m^^ + (25) 

+ (\A{\si-e)\^-e)-\A{\siTT-e)\^n-e)) . (26) 



This gives the Eq.(|TT]) for the quasiclassical energy E{Lz, 6) and the Eq. (|T^) 
for the shift Lf\e). The Eq.(|2§) is E(0, ^) = uo{e)Lf\e): it is the energy of 
quasiparticle moving along trajectory with 6 = 0. The nonzero value of this 
energy comes from the imaginary part of the order parameter Im A(|s|, —9) 
on the trajectory. 

4 Boundary conditions for quantization of bound 
states in the vortex core in chiral supercon- 
ductor. 

Let us consider simplest case of the axisymmetric vortex in p-wave supercon- 
ductor with the ^He-A order parameter. The order parameter in this vortex 
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is c{r){px + ipy)e^^'^, where c = ^ and Aq is the gap amphtude. The Bogoh- 
ubov Hamihonian for quasiparticles in the vortex core along the trajectory 
with given 9 is 

A (r) 

n = p-VT3. + Tiq^cir) cos(9 + N(f)) + T2q±c(r) sinW + Ncj)) , c(r) = °^ ' . 

Pf 
(27) 

The coordinate transformation which orients the axis x along the trajec- 
tory q^: 

— = -i—nds + Tic{r)coi^{N(P-iN-l)d))+T2c{r)siYi{N(P-{N-l)d)) . 
q± m 

(28) 

The second transformation is intended to delete the dependence on 9: 

^e'^N-i)r,ei2^ ^(29) 

^iiN-l)r,e/2^^-^(N-l)r,e/2 ^ {--T^ds + nc{r) COS N(j) + T2c{r) sin N(f)) .(30) 



The new Hamiltonian is the same as for the A^-quantum symmetric vortex 
in s-wave superconductor, but the Berry phase is different, now it is (A^ — 
l)r3^/2. This leads to different boundary condition for the wave function, 
which according to Eq.(^) is 

<iJ{e + 2n) = {-lf+^^{e) (31) 

This leads to the following quantization of levels in the core of the p-wave 
vortex with odd A^: 

E„ = -nuo , odd . (32) 

Thus for A^ = ±1 vortex there is a zero energy level. 

In general case the Cooper pair has an angular momentum projection 
m. Before we considered two examples: the s-wave superconductor, which 
belongs to the class m = 0; and the A-phase-like superconductors, which 
belongs to class m = ±1. In the case of general m, the Berry phase and 
boundary condition for the wave function of fermions in A^-quantum vortex 
are 

jSJ — ■rri 

OBerry = rs^^ , ^ (9 + 2n) = {-if-"^^ (9) (33) 

This leads to two classes of the fermionic spectrum in the symmetric core: 
En = nuo ii N — m is even and E^ = {n + l/2)ct;o if A^ — m is odd. 
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5 Effect of impurity scattering. 



Elastic scattering on impurity causes the transition between different trajec- 
tories. In the hmit of low energy of the quasiparticle the impact parameter of 
of the scattered particle tends to zero and becomes smaller than the distance 
from impurity to the center of the vortex. If the size of impurity is small 
then the scattering of the low-energy quasiparticle occurs only between the 
trajectories along the line between the vortex and impurities, i.e. between 
G = ^imp and 6* = TT — ^^i^p []TB[ . The " Josephson coupling" between these two 



states is 

T^i^p = 2\e'^m{7i - 0i„,p)^*(^in,p) + 2Ae-^^^*(7r - e,^^)m{e,^^) . (34) 
The Schrodinger equation is now 

(35) 

Let us choose the position of impurity at 7r/2. Then the relevant trajec- 
tories, where the scattering occurs, are at ±7r/2 and one has 

- + 2Ae-* (<) - * (-0 + 2Xe-^6 (» + f ) * (f ) = EmO) , 

(36) 

with boundary condition 

^ (-tt) = ±^ (tt) (37) 

Here the sign -|- is for the p-wave = 1 vortex and — is for s-wave = 1 
vortex. The solution is 

^ = A,e'^' , -1<0<1, (38) 

^ = A2e'^^ , -n <e <-^ , (39) 

^ = Age^^^ , ^<e <7l . (40) 

The conditions across two 5-function potentials and the boundary condi- 
tion Eq.(p7|) give 3 equations for the parameters A: 

- 2cuoe"'^ {Ai - A2) = Ae'V^ {Ai + A3) , (41) 
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• TVE 



-iujoe^ (A3 - Ai) = Ae-*^e"*^ {Ai + A2) , (42) 

e'^As = ±e"'^A2 . (43) 

Solution of these equations give the energy eigenvalues: 



ttE 2ujoX . 

cos = -r sm 7 , s — wave , ( 44 1 

CUq UJq + A"' 

sm = —cos 7 , p — wave . (45) 



In the s-wave case the Eq.(|4^) is similar to Eq.(2.10) of Ref. 0: the spectrum 
in the presence of impurity has the double period AE = 2ujo and consists of 
two equdistant sets of levels. These two sets transform to each other under 
symmetry transformation E —>■ —E, which is the "CPT" -symmetry of the 
system. 

For the p-wave case the Eg . (^Sf) also gives two sets of levels with the 
alternating shift. But the two sets are not mutually symmetric with respect 
to -E = 0. This contradicts to the "CPT" -symmetry of the system. The 
only way to reconcile the Eq. (|45|) with the symmetry is to fix the phase 
7 = 71/2 of the " Josephson coupling" . Then the energy levels are E^^ — n'^^o■, 
i.e. the same as without impurities. Thus the same "CPT" -symmetry, which 
is responsible for the eigenstate with E = 0, provides the rigidity of the 
spectrum. 



6 Conclusion 

The Berry phase in Eq. (|33|) is instrumental for the Bohr-Sommerfeld quan- 
tization of the energy levels in the vortex core. It chooses between the two 
possible quantizations consistent with the CPT-symmetry of states in super- 
conductors: En = nojQ and En = {n + 1/2)ujq. 

We found that the two spectra remain intact if the small perturbation is 
added, which violates the axial symmetry of the order parameter. We also 
considered the effect of single impurity on the spectrum of bound states. We 
found that if in a pure superconductor the spectrum is E^ = uujq (an example 
is = 1 vortex in chiral superconductor with m = 1) then the impurity 
does not change this spectrum. On the other hand if the initial spectrum 
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is En = {n + 1/2)070, the impurity splits it into two series according to the 
Larkin-Ovchinnikov prescription 0, This rigidity of the spectrum must be 
taken into account when the effect of randomness due to many impurities is 
considered with introduction of new level statistics for the fermionic spectrum 
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I am indebted to P. Wiegmann who noticed that there should be a Berry 
phase whcih is reponsible for quantization and to M. Feigel'man for fruitful 
discussions. 
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